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Abstract. We derive the formalism for steady state nonequilibrium dynamical mean-field 
theory in a real-time formalism along the Kadanoff-Baym contour. The resulting equations 
of motion are first transformed to Wigner coordinates (average and relative time), and then 
re-expressed in terms of differential operators. Finally, we perform a Fourier transform with 
respect to the relative time, and take the first-order limit in the electric field to produce the 
quantum Boltzmann equation for dynamical mean-field theory. We next discuss the structure of 
the equations and their solutions, describing how these equations reduce to the Drude result in 
the limit of a constant relaxation time. We also explicitly demonstrate the equivalence between 
the Kubo and nonequilibrium approaches to linear response. There are a number of interesting 
modifications of the conventional quantum Boltzmann equation that arise due to the underlying 
bandstructure of the lattice. 



1. Introduction 

Dynamical mean-field theory was introduced in 1989 by Metzner and Vollhardt J] as a particular 
limit where the spatial dimension approaches infinity. They showed that the self-energy is local 
in this case, so the many-body problem simplifies to that of an impurity in a time-dependent 
field plus a self-consistency relation [2]. Since then, nearly all of the classic problems of many- 
body physics have been solved (see [5] for a review). The generalization of the formalism to 
the nonequilibrium case is straightforward, but involves a number of technical complications. 
In this chapter, we describe the steady-state limit, where we start our system in equilibrium 
as to ~~ > —oo, then we turn on an external field at a later time t' (t' > to) and take the limit 
where t' Q — > — oo. We consider only a spatially uniform (but possibly time-dependent) electric 
field here, and ignore all magnetic-field effects. The Hamiltonian is written as the sum of two 
terms: (i) a time-independent Hamiltonian Ti which describes the field-independent piece, and 
(ii) the time-dependent piece Tif(t) arising from the field dependence. We adjust the electron 
density with a chemical potential, by subtracting the term (jlM from the Hamiltonian. The work 
presented here is a natural extension of the discussion in Mahan's review article [I] to the case 
of dynamical mean-field theory. 

The nonequilibrium formalism works with the so-called contour-ordered Green's function, 
which is defined for t and t' on the Kadanoff-Baym contour [2] (see Fig. ^) in terms of the 
electron creation (c r ) and annihilation (c r ) operators at lattice site r: 
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Figure 1. Kadanoff-Baym contour for the time variables. The time-domain cutoffs are 
symmetric at ±i max and we take the limit i max — ► oo. T/ie arrows indicate the direction to 
take along the contour for the time- ordering operation. 



with S C {TL') = exp[— i J c dt"Ti' I (t")] the evolution operator for the time-dependent piece of the 
Hamiltonian in the interaction representation, 0{t) = exp[i(Tt — fiN)t]0 exp[— i{TL — fiJ\f)t] 
the time-dependent operator for any operator O, Z = Trexp[— fi(TL — fiftf)], (3 is the inverse 
temperature, and T c is the time-ordering operator along the Kadanoff-Baym contour ("earlier" 
times ordered along the contour are moved to the right in the operator expressions). We work 
in units where fi = c = a = 1 (a is the lattice spacing) . 

The contour-ordered Green's function satisfies the Dyson equation (which can be thought of 
as the definition of the contour-ordered self-energy E c ) 

G c ry (t,t') = Gf^M + J dr" Jdt" Jtt'"G?y,{tX)^ c {t\t'")G c r ,,y{t'",t') (2) 

with each time integral being performed over the Kadanoff-Baym contour. The Green's function 
G c0 is the noniteracting contour-ordered Green's function in the presence of the field. Combining 
the nonequilibrium perturbation theory rules described by Langreth 6 with the perturbation 
theory expansion of Metzner [7], shows that the contour-ordered self-energy is local in infinite 
dimensions, and hence nonzero only when r = r'. This is why the self-energy has no spatial label 
in Eq. @. Since our field is spatially uniform, all Green's functions are translationally invariant, 
so transforming to momentum space diagonalizes the spatial part of the Dyson equation (note, 
we work in the Hamiltonian or temporal gauge with a spatially uniform vector potential before 
making the theory gauge- invariant). 

Keldysh [S] was the first to note that if we are interested in the steady-state limit, as described 
above, then the Kadanoff-Baym contour extends from — oo to oo, and we do not need to consider 
the matrix structure associated with the piece of the branch that extends along the imaginary 
axis; instead, we take it into account by incorporating the term exp[— (3(7i — [iN)]/ Z into the 
operator averages to represent the initial equilibrium distribution of the system for times in the 
distant past. Next, we can re-express the Dyson equation, which involves an integral over the 
Keldysh contour, into an integral over the real-time axis, but at the price of making the Green's 
function and self-energy into a 2 x 2 matrix with the corresponding Dyson equation [3]: 
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where matrix multiplication is understood amongst the 2x2 matrices and we have Fourier 
transformed from real space to momentum space (the momentum for the Green's functions is 
k). The three labels (R, K, and A) on the Green's functions and self-energies refer to the 
names retarded, Keldysh, and advanced, respectively; these objects are formed via specific 
linear combinations of the contour-ordered Green's function 8 . In order to describe these 
linear combinations, we must first define four Green's functions from natural choices of the time 
variables. The time-ordered Green's function is the contour-ordered Green's function with both 
t and t' on the upper real-time branch; the anti-time-ordered Green's function has both t and 
t' on the lower real-time branch; the lesser Green's function has t on the upper and t' on the 
lower real-time branches; and the greater Green's function has t on the lower and t' on the upper 
real-time branches. Then the retarded Green's function is one-half the time-ordered minus the 
anti-time-ordered plus the greater minus the lesser Green's functions. The advanced Green's 
function is one-half the time-ordered minus the anti-time-ordered minus the greater plus the 
lesser Green's functions. The Keldysh Green's function is the greater plus the lesser Green's 
functions. There is a conjugate Dyson equation as well, which looks similar to Eq. (J3J), but with 
the order of the G°T,G term switched to GXG°. 



2. Retarded Green's function 

In this section, we will develop the formalism for the retarded Green's function. To start, we 
examine the equation of motion (EOM) for the noninteracting retarded Green's function in the 
presence of a uniform electric field. The field is described with a time-dependent vector potential 
A(t) and a vanishing scalar potential in the Hamiltonian gauge, where E = —dtA(t). This is 
carried out with a Peierls substitution k — ► k — eA(i) to yield 



^ t + /i-e k _ eA(t) ]G£ (M') = 5(t-t'), 
^ + M-e k _eA(*')l Gk° (MO = S(t-t'), 



(4) 



with e k = — t* lim^oo J2i=i coskj/\/d; the hopping matrix element between nearest neighbors 
satisfies t = t* j2\fd in the infinite-dimensional limit [Q. The next step is to take each differential 
operator in the square brackets in Eq. © and operate each on G k (t, t') using the relevant Dyson 
equation. The equation that involves a derivative with respect to t uses the Dyson equation 
in Eq. (j3J), while the one with respect to t' uses the "conjugate" equation. We then take the 
sum and the difference of the two equations, and transform from the time variables t and t' 
to the Wigner coordinates jllj that involve the average and relative times: T = (t + t')/2 and 
t re i = t — t' . The Green's functions are now expressed in terms of the Wigner coordinates, for 
example, G^(T, t m \). Finally, we make the assumption that the retarded self-energy is a function 
of the relative time only (which will be argued to be true a posteriori for the Falicov-Kimball 
model |12j below). The two differential equations become 
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where we have assumed a constant uniform electric field, so A(t) = — Et; the electric field 
points along the diagonal (1, 1, 1, ...) of the infinite-dimensional hypercube. We need to perform 
some simple manipulations to get these equations into their final differential form. First we 
shift i — ► t + T — i re i/2 in each integral, and then we let t — ► t Te \ — i in the first integral 
of both equations. Finally, we use the fact that any infinitely differentiable function of time 
F(T) satisfies the Taylor series operator identity exp[idr]F(T) = F(T + 1) to produce the final 
time-dependent equations in the Hamiltonian gauge: 



l() t xel +P- ^ e k+eE(T+t rol /2) - 7jek+eE(T-t rel /2) 



G R (T, t re \) 



8(t Tel ) + J diZ R (t) cosh Qftr) G R (T,t iel - t), 

idr — £k+eE(T+t rol /2) + e k+eE(T-t rcl /2) G R (T, t re \) = 

-2 J diZ R (t) sinh Q#r) G£(r, t Icl - t). (6) 

The result in Eq. ©, is an exact equation of motion satisfied by the steady-state Green's 
function. It is complicated because it has many orders of the differential operators due to the 
cosh and sinh terms. It is also not yet in the form that we can take the linear limit needed 
for the quantum Boltzmann equation. The next step is to introduce a Fourier transform with 
respect to the t re \ coordinate: 

£ R (t re i) = 7^ J dcoe-^Z R (u;), G R (T,t vcl ) = ±- J duje~^G R (T,u). (7) 

When this is substituted in, we replace t terms by id u (this is accomplished by acting the 
differential operator on the exp[— io;t re i] term in the Fourier transform and then integrating by 
parts to move the differential operator from the exponential factor onto the self-energy factor), 
and we note that cosh(ia) = cos(a) and sinh(ia) = isin(a). Finally, we make the equations 
gauge invariant. For our choice of field, this process is straightforward, and we simply make the 
coordinate change k — > k — eET |13j . This requires us to modify the derivative with respect to 
T via 8t —> dx + eE- Vj^ [1]. Note that the gauge transformation operation is simply an average- 
time-dependent shift of the momentum, so the (local) self-energy is unchanged. Furthermore, 
the gauge-invariant Green's function G r (T,lo) = G R _ eET (T,Lo) turns out to be independent of 
T; we use the tilde symbol to denote gauge-invariant Green's functions. This last step follows 
from the fact that the self-energy has no T dependence, and the (steady-state) noninteracting 
gauge- invariant Green's function (for the simple cosine band), expressed as |141 115j 
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has no T dependence either. In Eq. l|H]). J n (z) is the Bessel function of the first kind with integer 
order n, E is the magnitude of each component of the electric field along the hypercube diagonal, 
and WBioch = eJ5 is the Bloch oscillation frequency; note that this choice for the electric field 
produces the natural Bloch frequency, even though the norm of the field grows like y/d. Since 
both G R0 and T, R are independent of T, it is a simple exercise to show that the gauge-invariant 
version of Eq. ©, found by making the gauge transformation explicitly on the equation, leads 



to a gauge-invariant retarded Green's function G R that is also independent of T. The final EOM 
becomes 



/ \ e ^^Ri , ^BlocliN e k/VR/ "^Bloch < 

(w + n)GQ(u) - —G£(uj + — g— ) - — 

= 1 + cos(^,eE • V k )E fl (o;0<5k K 
eE • V k G R (u) - e k G R (u + ^) + e k G R (u; 

= 2sin(^eE-V k )£^V)G k V) , (9) 

with e k = -i* lim^oc £f =1 sin(kj)/\/d. 

The equations in Eq. Q are exact for dynamical mean-field theory. One can easily verify 
that the gauge-invariant noninteracting Green's function in Eq. (JSJ satisfies the first equation 
of motion after using the identity aJ„+i(a) + aJ„_i(a) = 2nJ n (a). Furthermore, if the field 
vanishes, so E = 0, one can see that the first equation is satisfied by G R = l/[uj+fi — H R (ui) — e k ]. 
The second equation also holds trivially in equilibrium, because E = 0. We linearize Eq. Q in 
order to find the first nontrivial parts which are needed for the quantum Boltzmann equation. 
The first equation is the same as in equilibrium {E = 0), with the same solution [3], that is, 
the retarded Green's function is unchanged by the presence of the field to first order in E. The 
linearized form of the second equation is 
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G r (uj) = 0, (10) 



with v k = V k e k the velocity operator. It is a simple exercise to show that the equilibrium form 
of the Green's function also satisfies this equation [just note that G Req (u) depends on k only 
through e k , so that V k G Req (uj) = v k xdG Req (u;)/de k }. 

The final issue that needs to be resolved is that the retarded self-energy has no average time 
dependence. This is straightforward to show for the Falicov-Kimball model ^2], but requires too 
many steps to include the full derivation here. The basic idea is that the dynamical mean field for 
the impurity (which we map onto in the infinite-dimensional limit) is average time independent 
if the retarded Green's function and initial self-energy are both average time independent. 
Extracting the new self-energy from the impurity problem in the dynamical mean field does 
not introduce any average time dependence because all equations and fields are time-translation 
invariant for the steady-state, so the new self-energy is also average time independent. Since 
we can start the algorithm with T, R = 0, which has no average time dependence, the iterative 
algorithm to solve for the retarded self-energy and Green's function never generates any average 
time-dependence, so those functions will depend only on u. It is not clear at this point whether 
this holds for other models, because it is possible that the procedure to determine the retarded 
self-energy from the impurity problem in the dynamical mean field may introduce some average 
time dependence into the functions. Such a complication is straightforward to handle, but for 
simplicity, we have not considered it here. 



3. Keldysh Green's function 

The derivation for the Keldysh Green's function is similar, but it is complicated by the fact 
that the gauge-invariant Keldysh Green's function may still have average time dependence in 
the steady state. This possibility is related to the phenomenon of Bloch oscillations ^H]) where 
the steady state in the absence of scattering is an alternating current, oscillating at the Bloch 
frequency, which clearly is a state that has residual time dependence even in the steady state. 



This result arises because we are on a lattice and have restricted ourselves to a single band. In 
the continuum, the Keldysh Green's function is independent of the average time in the steady 
state. On the lattice, if we consider multiple bands, then if there is tunneling between bands, 
mediated by the field, then the Bloch oscillations may also be quenched and the system may 
evolve into an average time independent system. 

Since we have translational invariance, the spatial part of our problem diagonalizes when 
we work in momentum space, which we do here. We begin with the Dyson equation for the 
Keldysh Green's function, which follows from taking the Keldysh component of the full Dyson 
equation in Eq. ©. We can do the same with the "conjugate" equation as well. Finally, we 
can eliminate the Keldysh Green's function on the right-hand side by careful manipulation of 
the retarded or advanced Green's function Dyson equations. Hence, there are three forms for 
the Dyson equation of the Keldysh Green's function. We write each down schematically, with 
the notation that two time-dependent objects adjacent to one another imply a time integration 
over the real axis (not the Kadanoff-Baym contour): 
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(11) 



The noninteracting Keldysh Green's function satisfies a simple differential equation, because 
it has no discontinuity when t = t' like the advanced and retarded Green's functions. The 
equations of motion are 



[idt + n-ek+eEt] G£°(t,t') 
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(12) 



Now the derivation proceeds exactly like it did for the retarded Green's function. We operate 
each of the differential operators in brackets onto the Keldysh Green's function, and use the 
relevant Dyson equations, and the EOM of the noninteracting retarded, advanced, and Keldysh 
Green's functions to simplify the equations (the advanced Green's function is equal to the 
complex conjugate of the retarded Green's function with the time coordinates interchanged). 
Then we go to Wigner coordinates and change variables in the integrals. The only difference 
from the derivation for the retarded case is that the Keldysh self-energy generically depends on 
the average time. After some long algebra, we arrive at 



'■a j. 1 1 

ZC Vl + M - 2 e k+eE(T+t rcl /2) - 2 £ k+eE(T-t rcl /2) 



G)t(T, *rel) 



1 



dt 



e 2 



*G«(T,t Iel - i)Z R (t) + e-^G A (T,t vcl - t) 



e&G*(T,t iel -i)} [e-H^^^(T,t)] +e>G k f (T, t rcl - t)Z A (i)} 

\ldT - ek+eE(T+t rcl /2) + £k+eE(T-t rol /2) G k (T, t Ie \) = J dt 

{ e -I^G*(T,i rel -i)E*(i) + [e-I^G A (T,t iel -t)] [e^E^T.t) 
-p-G£(T,t rel -t)] [e-^^£*(T,t)l - e >G^(T,t rcl -t)S A (t)} 



(13) 



which are exact equations satisfied by the Keldysh Green's function. 

Next, we introduce the Fourier transform of the functions of i and i re i — i, and replace 
any remaining factors by derivatives with respect to the corresponding frequencies, and then 



integrating by parts as necessary to move the derivatives off of the exponential terms and onto 
the Green's functions or self-energies. Finally, we perform the gauge-invariance transformation 
by shifting the momentum vector with respect to the average time. We use the tilde to denote 
gauge-invariant Green's functions; the self-energies are unaffected by the transformation. The 
final exact equations of motion for the Keldysh Green's function are 

{lo + h)G^(T,uj) - ^e k G^(T,uj + ^u B ioch) ~ ^ e k G k( T ,u ~ T^Bioch) = 

+e~ i ^ eE - V *-^ T ) qR {uj ) z k (T, lo') + e-i 9 *' (^+eE- v k ) q k (r? u) {J) 
i&rG^(T,u)+ieE ■ V k Gg(T,cj) 
-ie k G£ (T, lo + ^Bioch) + ieuG*(T, lo - ^w B ioch) = 

ef^'^+ eE - v ^G^(r,a;)E fl ( W + e^ a "' eE ' Vk - a ^)G^(a;)S K (r,a;0 
-e-^' eE - v *- d ^G%(u)Z K (T,u') - e-^'^^G^T,^^'), 

(14) 

where, in both cases, we must take the limit lo' — > lo. In Eq. (j!4|) . we used the facts that the 
gauge-invariant retarded and advanced Green's functions and the local retarded and advanced 
self-energies are all independent of the average time T. 

Note that in equilibrium, with E = 0, the equilibrium forms for the Keldysh Green's function 
and self-energy 

G* eq (oj) = -2i[2/(w) - l]1mGg(u), X Keq (uj) = -2i[2/(w) - l]Im£*(u;), (15) 

do satisfy Eq. I|14jl. Since the EOM for the Keldysh Green's function is a homogeneous equation, 
the distribution-function factor [2/(cj) — 1] comes from the initial conditions (or equivalently 
from the bare Keldysh Green's function) and can be viewed as a boundary condition for the 
differential equations. 

The easiest way to see that the Keldysh Green's function and self-energy may generically have 
average time dependence, is to start with the gauge-invariant noninteracting Keldysh Green's 
function in a field (for the simple cosine band): 

Gg°(T, to) = 2vri[2/( ek _ eET - n) - 1] £ J n (^j 5 (u + /x - ^eEr?j . (16) 

The gauge-invariance transformation, which shifts the momentum wavevector, makes the Bessel- 
function and delta- function pieces independent of T, but the Fermi- factor piece [/(e) = 
l/{exp(/3e) + 1}] now picks up average time dependence. [One might be surprised by the fact 
that the noninteracting Keldysh Green's function has the Fermi-Dirac distribution that depends 
on ek in it, since the equilibrium Green's function described above had the factor /(lo) in it — this 
occurs because when the field vanishes, either form is correct since the noninteracting Keldysh 
Green's function has a 6(lo + \x — ek) term in it, but when a field is present, the noninteracting 
Keldysh Green's function has an infinite number of 5 functions coming from the Wannier-Stark 
ladder, and the initial equilibrium condition explicitly requires the Fermi-Dirac distribution as a 
function of ek — /U not lo.] This T-dependence is simple though: the band structure is a periodic 
function of momentum, which implies that the average time dependence must also be periodic 
in T; so the idea of a steady state that periodically repeats is possible [indeed, one only gets the 



correct Bloch oscillations for the noninteracting system if the form in Eq. (|16|) is used]. We need 
to show that both the Keldysh self-energy and the gauge-invariant Keldysh Green's functions 
are also periodic with the same period. If so, then the T-dependence can be treated in a Fourier 
series. It is straightforward to show from the full Dyson equation, that if the noninteracting 
Keldysh Green's function is periodic in T, and the Keldysh self-energy is also periodic in T, then 
the full Keldysh Green's function is periodic in T. Examining the impurity algorithm, shows that 
the Keldysh dynamical mean field is periodic in T. If the resulting impurity Green's function 
and self-energy are also periodic in time, then the impurity solver maintains the periodicity. 
This clearly holds for the Falicov-Kimball model ^2]) an d likely holds for many other models as 
well. Note that this analysis says the most general steady state must be periodic with the Bloch 
period. Such functions also include the constant function, since a constant is periodic with the 
Bloch period. 

Since the average time dependence is periodic, with a period 2n / 'wBiochi we can expand the 
Keldysh Green's function and self-energy in a Fourier series as 

n 

Y> K {T,uj) = ^e-™- hT £ K (n,w). (17) 

n 

Substituting these expansions into Eq. (|14|) produces essentially the same equations except the 
Keldysh functions depend on the Fourier index n and frequency uj, and the derivatives with 
respect to T are replaced by dr — > — inuiBioch- Finally, we take this equation and linearize it to 
obtain the final two expressions for the quantum Boltzman equation in dynamical mean field 
theory: 

(w + /x-e k )G , ^(n,tj) = G^(n,u;)ReS R (cj) + ReG R (oj)E K (n,oj) 
+ % -nuj moch G£(n 1 uj)d UJ lmZ R (uj) - ^eE • V k G^(n, u)dJmZ R {uo) 

+~nu mo ^dJmG R {u)Y. K {n,u) + ^eE • V k ImG R (u)d^ K (n, u), 

nwRiochGk (™> w ) + ieE • v kGk (n,u) - iek^Bioch^Gk (n, u) = 
2iGf(n,w)ImS fi (w) - 2ilmG R (u)E K (n, u) 
+nw B iochG^(n, w^ReE^u;) + ieE • V k G^(n, w)d w Re£ fl (w) 
-nuj B i och d u ReG R (uj)^ K (n,uj) + ieB ■ VkReG^i^d^ (n,u;). (18) 

These equations are obviously quite complex; the second equation is the one most closely related 
to the Blotzmann equation, and is the most important one for determining transport. 

As a first check, however, we evaluate these equations in equilibrium, with E = 0; in this case 
there is no average time dependence, so we are only interested in the n = equations. Taking 
the second equation, solving for G K , and substituting into the first equation for the G K term 
that appears on the right hand side, produces the relationship 

(u + fi — £k )G£ (c) = ( W + „ - eO [u + ^ _ _^ „ + [lm ^ R(u)]2 , (19) 

where the factor in the parenthesis cancels on both sides. After the cancellation, we are left 
with an expression that is consistent with the relation between the Keldysh self-energy and 
Green's function in equilibrium, but it does not explicitly produce the Fermi-Dirac distribution 
terms, as discussed above. The quantum Boltzmann equation typically is employed to obtain 
a relation between the self-energies and the Keldysh Green's function; it is expected that the 



self-energies will be determined by solving some kind of many-body theory that describes the 
system of interest (see the next section). 

Because the retarded and advanced Green's functions are unchanged to linear order in the 
field, and because we are only interested in the terms that are zeroth or first order in E for 
the quantum Boltzmann equation, we can replace the full Green's functions and self-energies 
by their equilibrium values when they are multiplied by either E or weioch [I]* This procedure 
produces a significant simplification of the results. For example, since the equilibrium Keldysh 
Green's function and self-energy are both independent of T, only the n = component survives, 
so all n^Bioch terms vanish. 

We perform the substitution of the equilibrium terms into their respective places, and then 
undergo some signficant algebra to end up at the final expression for the quantum Boltzmann 
equations: 

(u + fi- e k )Gg{n,uj) = Gg(n, w)Re£ Re<? (cj) + ReGf q {u;)Z K (n, u) 
-4ieE • vA/( W )[ImGf ? ( W )] 2 [ W + fl - ReE^V) - e k ]<5 n0 , 



[u + p- ReE^H - e k ] 2 + [ImE^u;)] 2 
-4ieE • v k [ImGf V)] 2 w /(w)<W (20) 

Note that the G£ and T1 K terms are not the equilibrium functions, but contain the equilibrium 
part plus the first-order shift in E; hence the first-order shift in the Keldysh Green's function 
arises from two sources — (i) the term coming from the electric-field dependence of the Keldysh 
self-energy and (ii) the explicit term given at the end of the second equation above (when n = 0) . 

The first equation is identical to the second multiplied by [to + \i — ReY, Req (uj) — e k ], so there 
is only one independent quantum Boltzmann equation. We use the second equation in Eq. (|2Uj) . 
because that is the equation most closely related to the Boltzmann equation (for a discussion of 
this issue, see Note further, that if the Keldysh self-energy has no average time dependence 
to first-order in E, then neither does the Keldysh Green's function. 



4. Linear-response conductivity in dynamical mean field theory 

The current operator is determined by the commutator of the polarization operator with the 
Hamiltonian, and can be expressed by 

J = e £ V k e k c k c k , 

k 

with c k (c k ) the electron creation (annihilation) operators for electrons with momentum k; this 
is the current per spin, since we are considering spinless electrons throughout this manuscript — 
adding spin increases the current by a trivial factor of 2. The term V k e k = v k is the electron 
velocity. We can calculate the expectation value of the current operator by employing the lesser 
Green's function G < = (G K — G R + G A )/2 as 

<J(T)> = 4ej^£ Vk G<(T,,,). (22) 

k 

Using the equilibrium forms for the retarded and advanced Green's functions (valid through 
first-order in E) and the quantum Boltzmann equation in Eq. (120H allows us to solve for the 
linear-response current 

<JCO> =e 2 £( E - v k)vk J ^[-^/(u,)][ImGfV)] 2 , (23) 



(21) 



and then the conductivity J2/3 a a^f3 = J« becomes 

a a p = 5 a ^e 2 ^l J ^[-^/M][ImGf V)] 2 , (24) 

k 

which is the dynamical mean field theory dc conductivity Note that in our derivation, 

we used the fact that the Keldysh self-energy is local and has no momentum dependence, so 
it vanishes when multiplied by the velocity and integrated over momentum, and only the term 
explicitly proportional to E in Eq. (|2U[) survives the summation over momentum. We did not 
need to assume that there was no average time dependence to the Keldysh Green's function, 
but if there is average time dependence, it arises from the Keldysh self-energy and cannot 
have any momentum dependence. This derivation is much simpler than that given in Mahan's 
review article 0] because the vertex corrections vanish in infinite dimensions ^Sj) so there is no 
Bethe-Salpeter equation to solve for the current; it simply follows directly from the quantum 
Boltzmann equation. Note that the full proof of the equivalence between the Kubo formula and 
the quantum Boltzmann equation was originally proved by Chen and Su who found an 
additional term neglected by Mahan, that is required for establishing the equivalence. 



5. Recovering the Drude-Sommerfeld model for transport 

The well-known expression for the Drude conductivity is: 

J ne 2 T 

&dc = = , (25) 

L m 

where n is the electron density, m is the effective electron mass, and r is the (constant) relaxation 
time, introduced phenomenologically by Drude. 

Here we concentrate on the conductivity of a strongly correlated system at finite temperature; 
all scattering is due to the electron-electron interaction. Determining the dc conductivity in this 
case is nontrivial (see, for example |20j). We analyze the behavior of the dc conductivity in the 
Falicov-Kimball model for spinless electrons on the infinite-dimensional hypercubic lattice. The 
Hamiltonian can be written in the following form |12j : 

H = ~lfn E c h + UJ2 W 4<H, (26) 
2Vd i 

where cj and c i are creation and annihilation operators for an itinerant electron at site i; t* = 1 
is the scaled nearest-neighbor hopping that we use for our energy unit; Wi is the localized 
electron number operator (equal to or 1); and U is the local Coulomb repulsion between 
localized and itinerant electrons. All electrons are spinless in this example. We also assume that 
the system is at half-filling, i.e., the particle density for the localized and itinerant electrons is 
equal to 1/2. 

In the limit of vanishing electron-electron interaction on the hypercubic lattice, an arbitrary 
constant electric field directed along the lattice diagonal (1,1,1,...), produces the current J(T) 
that can be calculated analytically from Eq. ()22j) . This current displays the expected Bloch 
oscillations [T5| : 

J(T) = -^-= sin^BiochT) J de^^p(e), (27) 

where p(e) = exp(— e 2 )/ v / 7r is the free-particle density of states, and the frequency of the 
oscillations is exactly the Bloch frequency weioch- It is interesting to discover whether these 
oscillations survive in the presence of electron-electron scattering, as U is turned on. It is 



difficult to observe this effect experimentally, since the period of oscillations for realistic fields 
is much larger than the typical scattering time, so the electron wavevector is randomized by 
scattering before it can undergo a Bloch oscillation. In the linear-response regime, the free- 
electron conductivity grows with the average time, since it is proportional to ET. The Bloch 
oscillations are absent because the period (proportional to 1/E) diverges as E — > 0. 

To calculate the expression for the dc conductivity from Eq. (|24[). one needs to find the 
equilibrium retarded Green's function for the interacting system. The case of the Falicov- 
Kimball model in the dynamical mean-field theory limit can be solved exactly I21| . This 
involves solving the set of the equations for the local Green's function G^ e c q (uj), the local self- 
energy Y, Req (uj), and the dynamical mean-field X Req (uj): 

G£V) = f dep{e) 1 — fj pn/ - ^ , (28) 



uj + p - e - Z Req (uj)' 
A* e » =u> + (j, — [G-^Vr 1 - £ Re9 M, (29) 



n Req, \ _ 1 ~ Wl Wl 

G loc ^> ~ u + ^ _ X Re q{uj) + u + ^ _ V _ X Re q{u) > 

where u?i is the average particle number for the localized electrons. In the case of half-filling we 
have Wi = 0.5 and \x = 0.5U. 

Once the retarded Green's function is found, the dc conductivity can be calculated from 
Eq. by introducing a joint density of states 



P2 



(e,e-) = _Lexp(-e 2 -e 2 ). (3 i) 



for the two energy functions and on the infinite-dimensional hypercubic lattice |22l I15j and 
integrating over e and e. 

The integration over e in Eq. (|24|) is straightforward and it gives the following general 
expression for the conductivity |17j : 



Vdc = o" / du}[-d u f(u)]r(u), (32) 

where o"o = e 2 ir/2d and the relaxation time satisfies 

/• [ImS Re 9(q;)] 2 /27r 2 

r(w) - y dep{e)j^ + ReJ , Req{oj) _ e]2 + [ ImS fe g(a;)]2} 2 • ( 33 ) 



It is possible to analyze the expression for adc analytically in the limit of low temperatures |23j . 
In this case the derivative —d u f{uj) is a function with a sharp maximum around lo = 0, and it is 
enough to approximate the relaxation time t(u>) by its low-frequency expansion. In the metallic 
case when U < U c = v2 it is sufficient to use just the lowest-order expansion t(u>) = t(0), 
yielding 

Ode = ^or(O) (34) 

which is the Drude form. The relaxation time is determined by the zero- frequency imaginary part 
of the retarded self-energy 

IniE^O), 

since the real part vanishes HeT, Req (u>) = (1 — Z)u> — > 
as uj — > (note that Z < for the Falicov-Kimball model because it is not a fermi liquid). In 
this case, the zero-frequency relaxation time can be evaluated as: 

r(0) ~ p(u) I dc [ ImSfleg (°)] 2 /^ 2 _ M 1 m) 



Therefore, the dc conductivity is inversely proportional to the imaginary part of the retarded 
self-energy in the metallic phase; this is the standard Drude picture. 

The situation is more complicated in the insulating case with U > y2. The interacting 
density of states Pi n t{w) = — (1 / ir)G^ e c q (uj) develops a pseudogap at u = 0. The real part of the 
self-energy is large for small since ReY, Req (uj) has a pole at u = 0. The imaginary part of 
the self-energy is small in this case except at u = where it diverges. In order to analyze the 
dependence of the dc conductivity on the model parameters in this case, we expand the local 
Green's function in powers of l/ReE^w) and ImE & "(w) in Eq. £8). The solution of the 
system of equations (|28 [1 -(|3U [1 for the self-energy to lowest order is [25] : 

ReS^H = (36) 
4a; 



u 2 



ivr[C/ 2 -2]<5H. (37) 



4a; 

The lowest-order approximation gives the following expression for the relaxation time |23| : 

16a; 4 

Substituting this expression into Eq. (|32j) yields (T = 1/(3 is the temperature here): 

T 4 4 f°° x 4 T A 7ir 2 

a ic = ^^j^-.Jjx—^ = vojfj^^. (39) 

The conductivity in the insulating phase goes to zero as T 4 , not exponentially, as expected for a 
conventional insulator with a true gap. This arises from the fact that there are a small number 
of states at low energy, but they have an exceedingly long lifetime, so they can carry a significant 
amount of current. 



6. Nonlinear response 

We have derived a series of exact differential equations that are satisfied by the nonequilibrium 
Green's functions in the steady state when a uniform constant electric field is applied [Eqs. © 
and ()13[)]. These equations appear to be intractable for a direct solution, because they involve 
differential operators to infinite order, and one needs to solve them with alternative strategies. 
Instead, we examined the linear-response regime, and found that we obtain the exact dc 
conductivity for dynamical mean field theory that was derived independently from the Kubo 
formula ^7j- One can ask, what happens when we go beyond the linear-response regime 
and examine nonlinear response? In the noninteracting case, the system will develop Bloch 
oscillations, so the presence of a dc field creates an ac current response. As the scattering is 
turned on, a classical Boltzmann equation analysis predicts that the oscillations will survive for 
a while, and then ultimately decay due to the scattering. Is it possible that in the steady state 
there are still oscillations, just with a reduced amplitude, or does the steady state always produce 
a constant current response? The answer to this question lies at the heart of whether or not the 
Keldysh Green's functions pick up average time dependence proportional to Vk (in a periodic 
fashion) for the steady state, or whether they have no average time dependence proportional to 
Vk (as in our linear-response example). There does not appear to be any simple way to resolve 
this question, but it seems possible to us that the Bloch oscillations can survive with a reduced 
amplitude, until the scattering gets to be too strong, and the system is forced into a constant, 
linear-response regime, due to the extremely short lifetimes. If this happens, it would be a 



truly quantum effect, since it does not seem to occur in the semiclassical Boltzmann equation. 
We intend to work further on this problem by directly investigating the solution of the Green's 
functions and self-energies in the steady state using techniques other than those derived from 
the equations of motion that have infinite-order differential operators in them. 

7. Conclusion 

In this work we have generalized the quantum Boltzmann equation approach of Mahan 0] to 
the case of dynamical mean field theory in a single-band model. This approach provides us 
with many simplifications, such as the fact that the self-energy has no momentum dependence, 
hence all momentum derivatives vanish. But it is complicated by some additional nonlinear 
effects brought on by the presence of the periodic band structure, which ultimately relate to 
the question of whether or not the system has Bloch oscillations, and how they evolve as the 
scattering is made stronger. Since we examined a linear-response limit to derive and use the 
quantum Boltzmann equation, we are not able to directly answer such questions here, instead 
we will do so in future publications. 

We were able to show that taking the linear limit of the nonequilibrium equations of motion, 
coupled with a simple theory for scattering by static charge defects, produced a quantum 
generalization of the semiclassical theory for transport with similar results to those of the Drude- 
Sommerfeld model. Along the way, we derived a series of exact differential equations that 
the nonequilibrium Green's functions satisfy, we showed how to incorporate periodic average 
time dependence into this framework if it exists in the system, and we established an exact 
equivalence between the quantum Boltzmann equation approach and the Kubo approach for the 
dc conductivity. 
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